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Abstract 

A method employing Bayesian statistics is used to incorporate recent experimental results on 
BjBj and BgBg mixing into a measurement of the Cabibbo-Kobayashi-Maskawa matrix elements 
with small theoretical uncertainties. The neutral B meson mixing results yield a slight improvement 
in the estimate of |Vtd|- Prospects for improving the knowledge of the CKM matrix elements with 
measurements of BgBg mixing and the K + — > tt + vv branching ratio are considered. 



The precise determination of the elements of the Cabibbo-Kobayashi-Maskawa (CKM) matrix jl| 
is one important goal of current experiments at BNL, CERN, CESR, FNAL, and SLAC and future 
experiments such as BaBar, Belle, CLEO-III, HERA-B, KTeV and LHC-B. Buras g § has pointed 
out that measurement of CP asymmetries in neutral B meson decays and the Kl — > ir vi> branching 
ratio can determine the CKM matrix with almost no theoretical uncertainties if V us is known. As 
yet, no measurements of these quantities exist. However, current experiments do provide information 
on two other useful measurements which are estimated to have 0(10%) theoretical uncertainties: the 
fractional B°B° mass difference and the K + — > tt + vv branching ratio. 



1 B°B° mixing 

In the Standard Model, mixing of neutral B mesons occurs via box diagrams dominated by internal 
top quark loops which allow the ratio of the CKM elements Vt s /V t d to be determined from B[|B[| and 
BgBg mixing measurements [||, 
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where Ama (Am s ) is the mass difference between the B[J (B|?) mass eigenstates, m s (rrtd) is the 
(Bg) mass and £ = 1.16±0.10 jjj is the ratio of hadronic matrix elements for the Bg and mesons and 
constitutes the theoretical uncertainty due to SU(3)fl avoU r breaking effects. Currently, B°B[] mixing 
is rather well measured, Am^ = 0.457 ± 0.019 ps -1 §; while experimental information from LEP on 
BgBg mixing indicates that Am s is at least ten times larger than Amd |], |7[ 10 1. Although the 
LEP experiments are incapable of resolving values of Am s greater than 10 ps -1 due to their limited 
event samples and proper time resolution, we show that it is nonetheless possible to determine the 
additional constraints placed on the CKM matrix elements from the experimental information on B 
meson mixing using Bayes' theorem. By way of demonstration, we first employ Bayes' theorem to 
determine the allowed range of Ara s given equation [l] and prior measurements of the magnitudes of 



CKM matrix elements [11]. 



For a proposition of interest q, prior knowledge i and additional information e, Bayes' theorem 

P(q\i A e) = P(q\i) P(e\i A q)/P(e\i) (2) 



allows one to adjust the "prior probability" of q, P(q\i), given the additional information e |Tz 
The factor P(e\i A q) is called the "likelihood" and often factorizes further since e is frequently the 
result of a sequence of statistically independent measurements and P(e\i) can often be determined by 
normalization [13, 14], Hfl . In our case, we use Bayes' theorem to demonstrate how probability densities 



related to the CKM matrix are affected by recent B°B° mixing results and how these probabilities 
might further be affected by potential measurements of the K + — > n + vi> branching ratio. 

Using the "standard" parameterization of the CKM matrix as recommended by the Particle Data 
Group (PDG) JIlJ], the CKM matrix is determined by four angles G = (#12, #23, #13> #13)- Assuming a 
uniform prior for Q within [0,tt/ 2 ) 3 x [0,2tt), the new probability density for is 

6 

P{@) = l[g(v i ,tii,a i ) (3) 

i=l 

where g(v,fx,a) is the Gaussian distribution and V\, V2, ■ ■ ■ , vq are the six quantities listed in Table [j] 



Matrix element 


Magnitude 


v ud 


0.9736 ± 0.0010 


v us 


0.2205 ± 0.0018 


V cd 


0.224 ± 0.016 


V cs 


1.01 ± 0.18 


v ub /v cb 


0.08 ± 0.02 


V cb 


0.041 ± 0.003 



Table 1: Experimental determination of the magnitudes of six of the CKM matrix 



elements as compiled by the PDG [11 



expressed as functions on [0,7r/2) 3 x [0, 2ir) with independently measured values /ij, fj,2, ■ ■ ■ , A*6 and 
corresponding variances o"f,cr|, . . . ,<7g. Once we have a probability density 0(0) then the density 
ip{x) for any real function /(0), 

Hx) = J^de <p(e) 17(6) < x] , (4) 
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can be evaluated numerically []T^, 16]. 

As an example of the method described above, consider the recent ALEPH study of the Bg Bg mass 
difference. As shown in Figure |], the experimental likelihood function £(Am s ) reflects the inability 
of the experiment to distinguish large values of Am s and only a lower limit on Am s is reported ||. 
However, assuming equation |l] one can use the measurements in Table [l] together with measurements 
of Amd II, the B° meson masses [11 1 and a theoretical estimate of £ Q as prior information in Bayes' 



theorem to improve upon the ALEPH likelihood. Assuming that these measurements and theoretical 
estimates are independent and normally distributed, the prior probability density of Am s can be 
constructed from equation || given the constraints in Table [l] and is shown in Figure |l[ There is a 95% 
probability that Am s > 5.8 ps _1 from the prior probability density alone. If the prior probability is 
combined with the ALEPH likelihood using Bayes' theorem, then the lower limit improves to 7.5 ps -1 
as shown in Figure [|. 




Figure 1: The dotted curve shows the prior probability of Am s determined from 
the measurements in Table [l]. The thin solid curve is C(Am s ) the likelihood of Am s 
from the ALEPH experiment || and the thick solid curve is the overall probability 
of Am s determined using Bayes' theorem. The normalization is arbitrary. 



The influence of the B°B° mixing measurements on the CKM matrix elements can be determined 
if equation |] is assumed. Given the ALEPH likelihood C(Am s ) and assuming that Amj, £, m s and 
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rrid are independent and Gaussian with means and variances as quoted above, then equation ||| can 
be used to determine the probability density for [Vt s /V t d| or the likelihood L(@) of the ALEPH 
results given the CKM angles. Using this likelihood and the PDG prior P(O), the density ipij for the 
magnitude of CKM matrix element Vjj is determined by 



AM « x- / de p (®) L ( ) riVv(e)| < «i 



(5) 



where the Vy(0) are the matrix elements expressed as a functions of the CKM angles. The resulting 
densities for the nine CKM elements, both with and without the B°B° mixing measurements, are 
shown in Figure ^. In the absence of B°B° mixing measurements, the 90% confidence limits on the 



matrix element magnitudes [17] are 



/ V u a V us V u b 
V c d V cs V c b 
^ V t d v ts V t b 



/0.9745 to 0.9757 0.219 to 0.224 0.002 to 0.005 
0.219 to 0.224 0.9736 to 0.9749 0.036 to 0.046 
v 0.006 to 0.013 0.035 to 0.045 0.9989 to 0.9993 



There is good agreement between these results and the PDG [[lTl for all elements except for Vt s , 
which is due to their rounding procedure [[18|], and V t d, which is obviously non-Gaussian. The PDG 



1 11, 15] and others [19, 2C] use the method of least squares to calculate confidence levels which can 
lead to inaccurate results if the predicted distribution of the magnitude of a matrix element is non- 
Gaussian. The two-lobed structure of | Vtd I is due to the phase 613 which is virtually unconstrained 
by the measurements of Table |l[ 




Figure 2: The probability density functions for the magnitudes of the nine CKM 
matrix elements. The solid curves show the probability densities determined with 
the PDG information in Table [1] alone; the dashed curves show the densities with 
the addition of the B°B° mixing results. 
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Only |Vtd| is significantly changed by the addition of the B°B° mixing information. As shown in 
Figure g, the 90% confidence limits change from (0.0057 to 0.0131) to (0.0054 to 0.0121) as large 
values of |V td | are disfavoured by £(Am s ) as shown in Figure |l[ 



2 The K + — > Ti + pp branching ratio 



A Bayesian estimate of the CKM matrix can also be constructed from K + — ► tt + ui> branching ratio 
measurements together with additional theoretical assumptions. The K + — » tt^v^^ branching ratio 
for a single lepton species I is given by pl| 



a 2 Br(K+ 



7T e 1 f, 



V* s V cd A^ L + V t *V td A(x t ) 



2 



Br(K+ - ir + v t v e ) = " . 4n x — - 2 , (6) 



2tt 2 sin 4 W | V u 

where xt = mf/mL. From Table 1 of reference [21|, we assume X^ NL = (11.2 rb 1.0) x 10 -4 for i = e,/i 
and X^ NL = (7.6 ± 1.0) x 10~ 4 for I = r where the quoted errors are dominated by uncertainties in 
Aj^g and the charm quark mass. We use the approximation X(xt) = 0.652^' 575 which is accurate 
to better than 0.5% for 150 < m t < 190 GeV/c 2 §], and take a(m w ) = 1/128, sin 2 W = 0.23, 



m t = 180 ± 12 GeV/c 2 [22||, m w = 80.32 ± 0.19 GeV/c 2 and Br(K+ -» vr°e + ^ c ) = 4.82 ± 0.06 % [1J. 



The probability density of i? = Br(K + — > 7t + z/j/), the branching ratio summed over £ = e, fx, r, can 
be determined using the procedure of equation || 

4>(B) = J q de P(9) r/?(9) < S] (7) 

where is now expanded to include xt, Br(K + — > TT°e + v c ) and X^ L as well as the six CKM mea- 
surements in Table [I]. We assume that xt, Br(K + — > 7r°e + f e ) and X^ L are independent and normally 
distributed about their central values but that the values of Xf^ L are completely correlated for the 
three lepton species. With these assumptions and the constraints of Table [l], there is a 90% probabil- 
ity that 0.5 x 10 -10 < Br(K + — > tx^vu) < 2.2 x 10~ 10 and the most probable value of the branching 
ratio is 0.8 x 10~ 10 . 

The current experimental upper limit at 90% confidence level of 2.4 x 10" 9 |23| from BNL E787 
is thus roughly an order of magnitude higher than the Standard Model prediction. From the rec- 
ommended procedure for the calculation of confidence levels for Poisson processes with background 
|nj, the probability density for Br(K+ -> ■k+vv) is P e (B) oc exp(-BAN K +) where A = 0.0027 is the 
experimental acceptance and A^ K + = 3.49 x 10 11 is the number of stopped K + [23]. The experimental 
result P e (B) is compared to the expected range from equation |7| in Figure ||. We also show the result 
of Bayes' theorem for the combined probability for the K + — * it^vv branching ratio. As expected 
the experimental result has no significant effect on the combined probability. The methods used in 
the paper should be useful in the near future as an upgraded version of BNL E787 is expected to 
be capable of observing the K + — ► ■n Jr vv decay assuming the Standard Model is correct [p^p . In the 
following section we describe the possible effect of such a measurement on the unitarity triangle. 



3 The unitarity triangle 

The unitarity triangle is a convenient way to present the relations between the CKM matrix elements 
\U, f§. The lengths of the sides of the triangle are given by (V* b V ud )/(V* b V cd ), (V* b V td )/(V* b V cd ) 
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Figure 3: The dotted curve shows the probability density of the K + — > 7r + //z/ 
branching ratio determined from the measurements in Table |1] and equations |6| and 
[7[ The thin solid curve is the probability of Br(K + — > ir + uu) determined from BNL 
E787. The thick solid curve is the combined probability of Br(K + — > -k + vt)). The 
normalization is arbitrary. 



and 1. The side of unit length conventionally has endpoints at (0,0) and (1,0) in the complex plane 
so that the apex of the triangle is given by (p, rf) . 

A straightforward generalization of equation || allows us to calculate the probability density of 
p and 77, ip(p,rj). Figure 0(a) shows the probability contours of ip(p,rj) given the measurements 
in Table |l|. The addition of the B°B° mixing information produces a slight modification of the 
contours as shown in Figure ^b). Finally we show the contours given hypothetical measurements of 
Ara s = 12.5 ± 0.5 ps" 1 and Br(K + -> tt + uv) = (1.0 ± 0.1) x 10~ 10 in Figure |(c). With these two 
measurements, there would be a 90% probability that 0.0077 < |V t d| < 0.0106 and the most probable 
value of |Vtd| would be 0.0088. As seen in Figure ||(c), it would be possible to confirm the prediction 
of CP violation (77 7^ 0) by the Standard Model at the 95% confidence level with these two particular 
measurements. 

4 Conclusions 

Using Bayes' theorem, recent measurements of B°B° mixing are used to improve upon current esti- 
mates of the CKM matrix elements with minimal theoretical uncertainty. Only |V t( j| is significantly 
affected by these results. The impact of current and potential measurements of B°B° mixing and the 
K + — > ir + vi> branching ratio on the individual CKM matrix elements and the CKM unitarity triangle 
is also examined. 

This work was supported by US Department of Energy contracts DE-FG05-92ER40742 and DE- 
FC05-85ER250000. 
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Figure 4: The 68.3% (dark) and 95.5% (light) probability contours for the apex 
of the unitarity triangle determined using (a) the measurements listed in Table |I| 
only, (b) the B°B° mixing results and Table [1], and (c) Table [I] and hypothetical 
measurements of Am s = 12.5 ±0.5 ps _1 and Br(K + — > = (1.0 ±0.1) x 10 -10 . 

A non-zero r\ would confirm CP violation as predicted by the Standard Model. 
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